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We study the nonequilibrium phase transition in a model of aggregation of masses allowing for 
diffusion, aggregation on contact and fragmentation. The model undergoes a dynamical phase transition 
in all dimensions. The steady state mass distribution decays exponentially for large mass in one phase. 
On the contrary, in the other phase it has a power law tail and in addition an infinite aggregate. The 
model is solved exactly within a mean field approximation which keeps track of the distribution of 
masses. In one dimension, by mapping to an equivalent lattice gas model, exact steady states are 
obtained in two extreme limits of the parameter space. Critical exponents and the phase diagram are 
obtained numerically in one dimension. We also study the time dependent fiuctuations in an equivalent 
interface model in (1+1) dimension and compute the roughness exponent x and the dynamical exponent 
z analytically in some limits and numerically otherwise. Two new fixed points of interface fiuctuations 
in (1 + 1) dimension are identified. We also generalize our model to include arbitrary fragmentation 
kernels and solve the steady states exactly for some special choices of these kernels via mappings to 
other solvable models of statistical mechanics. 



I. INTRODUCTION 



By now, there is a fairly good understanding of the 
nature and properties of phase transitions in systems in 
thermal equilibrium, as one changes the strengths of ex- 
ternal fields such as temperature, pressure or magnetic 
field. On the other hand there is a wide variety of in- 
herently nonequilibrium systems in nature whose steady 
states are not described by the equilibrium Gibbs dis- 
tribution, but are instead determined by the underlying 
microscopic dynamical processes. The steady states of 
such systems may undergo nonequilibrium phase transi- 
tions as one changes the rates of the underlying dynami- 
cal processes. As compared to their equilibrium counter- 
parts, these nonequilibrium steady states and the tran- 
sitions between them are much less understood owing to 
the lack of a general framework. It is therefore necessary 
to study simple models of nonequilibrium processes, both 
in order to discover new types of transitions as well as to 
understand the mechanisms which give rise to them. 

In this paper, we study nonequilibrium phase transi- 
tions in an important class of systems which involve the 
microscopic processes of diffusion, aggregation upon con- 
tact and fragmentation of masses. These processes arise 
in a variety of physical settings, for example, in the for- 
mation of colloidal suspensions ^ , polymer gels , river 
networks aerosols and clouds They also enter in 
an important way in surface growth phenomena involv- 
ing island formation j^]. Below we introduce a simple 
lattice model incorporating these microscopic processes 
and study the nonequilibrium steady states and the tran- 
sitions between them both analytically within mean field 
theory and numerically in one dimension. Some of the re- 
sults of this paper have been reported earlier in a shorter 



version 

The paper is organized as follows. In section II, we de- 
fine the model and summarize our main results. In sec- 
tion III, we solve the mean field theory exactly and char- 
acterize the phases and the transitions between them. In 
section IV, we report the results of numerical simulations 
in one dimension for both symmetric and asymmetric 
transport of masses. In section V, the model in one di- 
mension is mapped exactly to a lattice gas model whose 
properties are used to deduce the steady state mass dis- 
tribution exactly in two extreme limits. In section VI, 
we map this lattice gas model further to an interface 
model and make connections to other well studied inter- 
face models in some limits. We study the dynamics via 
computing the width of the interface both analytically in 
some limiting cases and numerically otherwise. We iden- 
tify two new fixed points of interface dynamics in (1-1-1) 
dimension. In section VII, we generalize our model to 
include an arbitrary fragmentation kernel and obtain ex- 
act results for special choices of this kernel via mappings 
to other solvable models of statistical mechanics. In the 
Appendix we outline the exact solution for uniform frag- 
mentation kernels. We finally conclude with a summary 
and a few open questions in section VIII. 



II. THE MODEL 

Our model of diffusion, aggregation and dissociation is 
defined on a lattice, and evolves in continuous time. For 
simplicity we define it on a one-dimensional lattice with 
periodic boundary conditions although generalizations to 
higher dimensions are quite straightforward. Beginning 
with a state in which the masses are placed randomly. 
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a site i is chosen at random. Then one of the following 

events can occur: 

1. Diffusion and Aggregation: With rate pi, the mass 
rrii at site i moves either to site i — 1 or to site i + 
If it moves to a site which already has some parti- 
cles, then the total mass just adds up; thus ^ 
and mi±i — > mi±i + rrii. 

2. Chipping (single-particle dissociation): With rate 

P2, a bit of the mass at the site "chips" off, i.e. 
provided > 1, a single particle leaves site i and 
moves with equal probability to one of the neigh- 
bouring sites i — 1 and thus mj — > mj — 1 and 
mi±i — > mi±i + 1. 

We rescale the time, t p\t, so that the diffusion and ag- 
gregation move occurs with rate 1 while chipping occurs 
with rate w = p2/pi- 

This model clearly is a very simplified attempt to de- 
scribe systems with aggregation and dissociation occur- 
ing in nature. For example, if one is thinking of gelation 
phenomena, then a polymer of size k is represented by a 
point particle of mass k in our model. Thus we ignore 
the spatial shape of the real polymer which however can 
play an important role under certain situations. We have 
also assumed that the fusion of masses after hopping or 
chipping occurs instantaneously, i.e., the reaction time 
scale is much smaller than the diffusion time scale. Thus 
our model is diffusion- limited. A somewhat more severe 
assumption is however that both desorption and diffusion 
rates arc independent of the mass. In a more realistic sit- 
uation these rates will depend upon the mass. However, 
our aim here is not to study any specific system in full 
generality, but rather to identify the mechanism of a dy- 
namical phase transition, if any, in the simplest possible 
scenario involving these microscopic processes. If one is 
interested in a more realistic description of any specific 
system, one could and should include these features in the 
model. But for the purpose of this paper, we stay with 
the simplest version and show below that even within 
this simplest scenario novel, dynamical phase transitions 
occur which are nontrivial yet amenable to analysis. 

In this model, the total mass M is conserved and fixed 
by the initial condition. Let p = M/N denote the den- 
sity, i.e., mass per site where N is number of sites of the 
lattice. In the above definition of the model, mass at 
each site can move symmetrically either to the left or to 
the right with equal probability. We call this the symmet- 
ric conserved-mass aggregation model (SCMAM). In this 
paper we also study the fully asymmetric version of the 
model where masses are constrained to move only in one 
direction (say to the left). We call this the asymmetric 
conserved-mass aggregation model (ASCMAM). In the 
sections that follow, we show that the nonequilibrium 
critical behaviour of SCMAM and ASCMAM belong to 
different universality classes. This can be traced to the 
fact that in the asymmetric case there is a nonzero mass 
current density in the system. 



In both of these models, there are only two parameters, 
namely the conserved density p and the ratio w = p2 /pi 
of the rate of chipping of unit mass to that of hopping of 
the entire mass on a site, as a whole. The question that 
we mainly address in this paper is: given (p, w), does the 
system reach a steady state in the; long time limit? If 
so, how can we characterize this steady state? We show 
below that indeed for all (p,w), the system does reach 
a steady state. This is a nonequilibrium steady state in 
the sense that for generic values of (p, w) it is not de- 
scribed by the Gibbs distribution associated with some 
Hamiltonian. In order to characterize the steady state we 
study the single site mass distribution function P(m.,t) 
as the time i — > oo. We show below that there exists a 
critical curve in the (p, w) plane across which the steady 
state behaviour of the system, as characterized by P{m), 
undergoes a novel phase transition. 

Let us siimmarize our main results: In our model, there 
are two competing dynamical processes. The diffusion 
cum coalescence move tends to produce massive aggre- 
gates at the expense of smaller masses, and in this pro- 
cess, also creates more vacant sites. The chipping of sin- 
gle units of mass, on the other hand, leads to a replenish- 
ment of the lower end of the mass spectrum. The result 
of this competition is that two types of steady states are 
possible, and there is a dynamical phase transition be- 
tween the two across a critical line Pc{w) in the {p,w) 
plane. For a fixed w, ii p < Pc{w), the steady state mass 
distribution P{m) decays exponentially for large m. At 
p = Pciw), P{m) decays as a power law P(m) ~ for 
large m, where the exponent r is the same everywhere 
on the critical line Pc{w). A more striking and inter- 
esting behaviour occurs for p > pc{w). In this phase, 
P(m) decays as the same power law ^ m~'^ for large m 
as at the critical point, but in addition develops a delta 
function peak at m = oo. Physically this means that an 
infinite aggregate forms that subsumes a finite fraction 
of the total mass, and coexists with smaller finite clus- 
ters whose mass distribution has a power law tail. In the 
language of sol-gel transitions, the infinite aggregate is 
like the gel while the smaller clusters form the sol. How- 
ever, as opposed to models of irreversible gelation where 
the sol disappears in the steady state, in our model the 
gel coexists with the sol, which has a power law distribu- 
tion of mass. Interestingly, the mechanism of the forma- 
tion of the infinite aggregate in the steady state resem- 
bles Bose-Einstein condensation (EEC), though the con- 
densate (the infinite aggregate here) forms in real space 
rather than momentum space as in conventional EEC. 

This nontrivial dynamical phase transition occurs in 
both the models SCMAM and ASCMAM in all spatial 
dimensions d including d = I. We expect that the ex- 
ponent T depends on the dimension d but is universal 
with respect to lattice structures and initial conditions. 
However, the bias in the movement of masses that distin- 
guishes the two models SCMAM and ASCMAM is a rel- 
evant perturbation and the corresponding exponents Tg 
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(for SCMAM) and (for ASCMAM) differ from each 
other. 

We comment on the relationship of our model and re- 
sults to earlier work on related models. 

(i) Takayasu and coworkers have studied a lattice 
model where masses diffuse and aggregate upon contact 
as in our model. However our model differs from the 
Takayasu model in the following important way. In that 
model there is a nonzero rate of injection of a single unit 
of mass from the outside into each lattice site, whereas 
in our model the Takayasu injection move is replaced by 
the 'chipping' of a single unit of mass to a neighbouring 
lattice site. Thus in our model total mass is conserved 
as opposed to the Takayasu model where total mass in- 
creases linearly with time. In the Takayasu model the 
mass distribution P(m) has a power law decay in the 
steady state but there is no phase transition as in our 
conserved model. Finally, while the directed and undi- 
rected versions of the Takayasu model evolve in the same 
way, in our model directionality in motion changes the 
universality class. 

(ii) In the context of polymer chain growth kinetics, mod- 
els of aggregation in 'dry' environments have been stud- 
ied earlier only within a rate equation approach [^,^. 
Although the aggregation of polymers and dissociation 
of single monomers were allowed in those models, they 
lacked the important process of the local diffusion of 
masses which is included in our lattice model. 

(iii) Models of vacancy cluster formation consider the at- 
tachment and detachment of single vacancies to clusters 
I pof , often with rates that depend on the cluster size. 
This would correspond to allowing only chipping moves 
in our model (Section V A). Since cluster aggregation 
moves are absent in this case, there is no tendency to 
form very large clusters and the mass distribution decays 
exponentially with m. 

(iv) Recently, lattice gas models have been proposed ||Tl| ] 
to describe the distribution of droplets in fast-expanding 
systems such as in the fragmentation process following a 
nuclear collision. The distribution of fragments shows a 
pronounced peak at the large mass end, reminiscent of 
our infinite aggregate. However, the distribution of the 
remaining fragments decays exponentially, and not as a 
power law as in our case. 

(v) Bose-Einstein-like condensation in real space has also 
been found in lattice gas models of traffic in which 
different cars have different maximum speeds, chosen 
from an a priori specified distribution. This corresponds 
to phase separation into high density and low density 
regions, a phenomenon which is also found when ran- 
domness in hopping rates is associated with fixed points 
in space, rather than with cars [ p^Jl5| . An important 
difference between these studies and ours is that they 
involve quenched disorder (the assignment of different 
maximum speeds, either to cars, or to different lattice 
sites), whereas the BEC phenomenon in our case occurs 
in the absence of disorder, in a translationally invariant 
system. 



III. MEAN FIELD THEORY 



In this section, we study the conserved mass models 
within the mean field approximation which keeps track 
only of the distributions of masses, ignoring correlations 
in the occupancy of adjacent sites. The mean field the- 
ory is identical for both the symmetric and asymmetric 
models. This is a defect of the MF approximation as, 
in fact, the existence of a nonzero mass current in the 
asymmetric model affects fluctuations in an important 
way. Although MF theory misses this important aspect 
of the physics, it is still instructive since it reproduces the 
phase diagram correctly, at least qualitatively. Besides, 
in high dimensions where fluctuations are negligible, the 
MF theory captures the correct physics even quantita- 
tively. 

Ignoring correlations between masses at neighbouring 
sites on the lattice, one can directly write down the evo- 
lution equation for P(m,t), the probability that any site 
has a mass m at time t. 



dP{m,t) 
Jt 

dP{0,t) 
Jt 



-(1 + + sit)]P{m,t) +wP{m + l,t) 



ws{t)P{m-l,t) + P*P- m>l 



(1) 



-(1 + w)s{t)P{0, t) + wP{l, t) + s[t). (2) 



Here s{t) = 1 — P{Q, t) is the probability that a site is oc- 
cupied by a mass and P*P — X]m'=i ^("^'7 t)P{m—m', t) 
is a convolution term that describes the coalescence of 
two masses. 

The above equations enumerate all possible ways in 
which the mass at a site might change. The first term 
in Eq. (|l|) is the "loss" term that accounts for the prob- 
ability that a mass m might move as a whole or chip 
off to either of the neighbouring sites, or that a mass 
from a neighbouring site might move or chip off to the 
site in consideration. The probability of occupation of 
the neighbouring site, s{t) = X]m=i ("^i multiplies 
P{m, t) within the mean-field approximation where one 
neglects the spatial correlations in the occupation proba- 
bilities of neighbouring sites. The remaining three terms 
in Eq. (|l]) arc the "gain" terms enumerating the num- 
ber of ways that a site with mass m' ^ m can gain the 
deficit mass m — m' . The second equation Eq. (H) is a 
similar enumeration of the possibilities for loss and gain 
of empty sites. Evidently, the MF equations conserve the 
total mass. 

To solve the equations, we compute the generating 
function, Q{z,t) X]m=i ^(™' *)-^™ from Eq. (|l|) and 
set dQ/dt = in the steady state. We also need to use 
Eq. (||) to write P{^,t) in terms of s{t). This gives us a 
quadratic equation for Q in the steady state. Choosing 
the root that corresponds to Q{z = 0) = 0, we find 
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w + 2s + ws w wsz 



+ ^^-^^ ^{z- zi){z- Z2). (3) 

where Zi^2 = (w + 2 =F 2y/w + l)/ws. The value of the 
occupation probability s is fixed by mass conservation 
which implies that ^ mP(m) — M / L = p. Putting 
dzQ{z = 1) = p, the resulting relation between p and 
s is 

2p = w{l - s) - ws^{zi - 1)(Z2 - 1) . (4) 

The steady state probability distribution P{m) is the 
coefficient of in Q{z) and can be obtained from Q{z) 
in Eq. (||) by evaluating the integral 

P{m) I ^dz (5) 

over the contour Co encircling the origin. The singular- 
ities of the integrand govern the asymptotic behaviour 
of P{m) for large m. Clearly the integrand has branch 
cuts at z = Z1.2. For fixed w, if one increases the density 
p, the occupation probability s also increases as evident 
from Eq. (|4|). As a result, both the roots zi_2 start de- 
creasing. As long as the lower root zi is greater than 1, 
Eq. (P is well defined and the analysis of the contour 
integration around the branch cut z — zi, yields for large 
m, 

P(m) - 6-"/"*/"^^^^ , (6) 

where the characteristic mass, m* = l/log(zi). It di- 
verges as ~ (sc — s)~^ as s approaches Sc = {w + 2 — 
2\/w -\- ^)/w. Sc is the critical value of s at which zi = 1. 
This exponentially decaying mass distribution is the sig- 
nature of the "disordered" phase which occurs for s < Sc 
or equivalently from Eq. (Q) for p < Pc{w) = \/w + I — 1. 

When p = pc, wc have zi = 1, and analysis of the 
contour around z = zi = 1 yields a power law decay of 
P(m), 

P(m) - (7) 

As p is increased further beyond Pc{w), s cannot increase 
any more because if it does so, the root zi would be less 
than 1 (while the other root Z2 is still bigger than I) and 
the right hand side of Eq. (^) would become complex. 
The only possibility is that s sticks to its critical value Sc 
or equivalently the lower root zi sticks to 1. Physically 
this implies that adding more particles does not change 
the occupation probability of sites. This can happen only 
if all the additional particles (as p is increased) aggregate 
on a vanishing fraction of sites, thus not contributing to 
the occupation of the others. Hence in this "infinite- 
aggregate" phase P(m) has an infinite-mass aggregate, 
in addition to a cluster distribution with a power law de- 
cay m^^/^. Concomitantly Eq. (^ ceases to hold, and 
the relation now becomes 



P = ^(1 - Sc) +P00 (8) 

where poo is the fraction of the mass in the infinite aggre- 
gate. The mechanism of the formation of the aggregate is 
reminiscent of Bose Einstein condensation. In that case, 
for temperatures in which a macroscopic condensate ex- 
ists, particles added to the system do not contribute to 
the occupation of the excited states; they only add to the 
condensate, as they do to the infinite aggregate here. 

Thus the MF phase diagram (see inset of Fig. 1) of 
the system consists of two phases, "Exponential" and 
"Aggregate", which are separated by the phase bound- 
ary, pc{w) = y/w + 1 — 1. We note that in conserved- 
aggregation models in 'dry' environment studied earlier 
within a rate equation approach l^],^, the steady state 
mass distribution also changed from an exponential dis- 
tribution to a power law as the density was increased 
to a critical value. However, the existence of the infi- 
nite aggregate in the steady state for p > Pc{w) was not 
identified earlier. 



IV. NUMERICAL SIMULATION IN ONE 
DIMENSION 

In order to see if the MF phase diagram remains at 
least qualitatively correct in lower dimensions, we have 
studied both SCMAM and ASCMAM using Monte Carlo 
simulations on a one-dimensional lattice with periodic 
boundary conditions. Although we present results here 
for a relatively small size lattice, N = 1024, we have 
checked our results for larger sizes as well. We confirmed 
that all the qualitative predictions of the mean-field the- 
ory remain true in 1-d though the exponents change from 
their MF values. 

Figure 1 displays two numerically obtained plots of 
P{m) in the steady state of SCMAM. For fixed w = 1.0, 
we have measured P{m) for two values of the density, 
namely p = 0.2 and p = 3.0. For p = 0.2, we find ex- 
ponential decay of P(m) (denoted by x in Fig. 1). For 
p — 3.0, P{m) decays as a power law (denoted by + 
in Fig. 1) which is cut off by finite size effects but in 
addition, there is a sharp peak at a much larger mass 
signalling the existence of the 'infinite' aggregate as pre- 
dicted in the MF theory. We confirmed that the mass 
Magg in this aggregate grows linearly with the size, and 
that the spread SMagg grows sublincarly, implying that 
the ratio SAIagg/Magg approaches zero in the thermody- 
namic limit. As one decreases p for fixed w, the mass 
Magg decreases and finally vanishes at the critical point, 
Pc(w) where P(m) only has a power law tail. For w; = 1, 
we find numerically Pc(l) ~ 0.39. In the inset of Fig. 1, 
we plot the numerical phase boundary in the (p, w) plane 
(denoted by closed circles) . This is only a rough estimate 
of the phase boundary obtained from small lattice sizes 
and the numerical accuracy of these points should not 
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be taken too seriously. For comparision, we also plot the 
MF phase boundary Pc{w) = {w + 1) — 1. 

In Fig. 2 and in its inset, we present similar plots for 
the asymmetric model ASCMAM. Here the steady state 
mass distribution function P{m) is plotted for two values 
of the density p = 0.2 and p = 10.0 with fixed w — 1.0. 
In the first case, P[m) decays exponentially whereas in 
the second case it has a power law tail and in addition 
the 'infinite' aggregate as in the case of SCMAM. 

The exponent for SCMAM which characterizes the 
finite-mass fragment power law decay for p > pdw) 
is numerically found to be 2.33 ± .02 and remains the 
same at the critical point p = pdw). In the asymmet- 
ric model ASCMAM, we find the corresponding expo- 
nent Tas ~ 2.05 within numerical error. Note however 
that because the total mass and hence the mass density, 
p — J2 mP{m) is conserved and finite, the decay of P(m) 
must be faster than for large m to avoid ultraviolet 
divergence. In ASCMAM the numerical value of Tas in 
1-d is very close to 2 suggesting perhaps that P{m) de- 
cays as with additional logarithmic corrections such 
that p remains finite. But within our simulations, it is 
not easy to detect these additional logarithmic factors. 
Thus clearly in 1-d, SCMAM and ASCMAM belong to 
different universality classes. 



V. MAPPING TO A LATTICE GAS MODEL IN 
ONE DIMENSION 

In this section, we show that in one dimension the mass 
model studied above can be mapped exactly onto a lat- 
tice gas (LG) model consisting of particles and holes. In 
the language of the LG model, it is somewhat easier to 
understand the two phases and the transition between 
them. Besides, for certain limiting values of the param- 
eters, the steady state of the LG model can be solved 
exactly. This then provides, via the mapping, exact so- 
lutions for the mass model in those limiting cases. 

Consider the mass model (both SCMAM and ASC- 
MAM) on a ring R oi N lattice sites. Let rrii denote 
the mass at site j of in a given configuration. Let 
M — X^ili denote the total mass on R. We first con- 
struct a new ring R' consisting oi L = N + M sites. For 
every lattice site i of ring R, we put a particle (labelled 
by i) on ring R' such that the i-th and (i -I- l)-th parti- 
cle on R' are separated exactly by rrii holes. The ring R' 
will therefore have N particles and M holes. Also by con- 
struction, these particles on R' are hard core, i.e., any site 
of R' can contain at the most one particle. Thus every 
mass configuration on R maps onto a unique particle- 
hole configuration on R' . In Fig. 3, we give an example 
of this mapping. We also note that particle density p' 
on R' is simply related to the mass density p on R via, 
p' ^N/{N + M) ^ 1/(1 + p). 

Given this exact mapping of a mass configuration on 
i? to a particle hole configuration on i?', we now exam- 



ine the correspondence between the mass dynamics on 
R and the particle dynamics on R' . Consider a pair of 
neighbouring sites (i — 1) and i on R where the masses 
are m^-i and rrii. This translates to having mi_i holes 
to the right of (i — l)-th particle and holes to the right 
of i-th particle on R' . Consider first the chipping move 
that occurs with rate p2 (Fig. 3). If a single unit of mass 
chips off the (i — l)-th site and moves to i-th site (i.e., 
rrii^i — > TOi_i — 1 and mi mi + 1) on R, it corresponds 
to the i-th particle on R' hopping to its neighbouring site 
to the left with rate p2. Similarly, the reverse move of 
a unit mass chipping off i-th site to (i — l)-th site on R 
(i.e., mi ^ mi — 1 and mi-i — > mi-i + 1) corresponds 
to the i-th particle on R' hopping one step to the right 
with rate p2 ■ Thus the chipping move of SCMAM on R 
corresponds precisely to the "symmetric exclusion pro- 
cess" (SEP) [|l^ on R' where a particle can hop to its 
nearest neighbour on either side provided the neighbour 
is empty. Similarly the chipping move of ASCMAM on 
R corresponds exactly to the "asymmetric exclusion pro- 
cess (ASEP) 11,0,111 on R' where particles move only 
along one direction on R' . 

But in addition to chipping, we also have the diffu- 
sion and aggregation move in the mass model that oc- 
curs with rate pi (as shown for the block of particles 
in Fg. 3). Consider once again a pair of sites (i — 1) 
and i with respective masses mi^i and TOj on R. Sup- 
pose the whole mass mi-i moves to the i-th site (i.e., 
mi_i — > and mi — > mi + m^-i). This would mean 
that on R' , the i-th particle jumps to the farthest avail- 
able hole (without crossing the (i — l)-th particle) to its 
left with rate pi. Similarly the reverse move, — > 
and mi_i —>■ mi^i + mi, would translate on R' to the 
i-th particle jumping to the farthest available hole to its 
right (without crossing the (i -I- l)-th particle) with rate 
Pi. In the asymmetric version of the model, the particle 
can jump to the farthest available hole (without crossing 
the next particle) in one direction only. 

To summarize, in our LG model we have a bunch of 
hard core particles with particle density p' . The world 
lines of particles can not touch or cross each other due 
to their hard core nature. There are two possible moves 
for each particle. With rate p2, a particle moves to its 
adjacent site (if it is empty) and with rate pi the particle 
jumps to the farthest available hole maintaining the hard 
core constraint. In SCMAM, the particle can hop with 
equal probability to the left or right. In ASCMAM, it 
hops only in one direction, say to the right. In this LG 
language, it is the competition between the short range 
and long range hopping of the particles that is responsi- 
ble for the phase transition. As in the mass model, the 
only two parameters of the model are the ratio of the two 
rates, w = P2/P1 and the density of particles p'. We also 
note that P{m) in the mass model corresponds to the 
size distribution of hole clusters in the LG model. 

In the following we fix the density of particles p' and 
study the two extreme limits w ^ 00 and w = Q where 
exact results can be obtained. 
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A. Only chipping : ui oo 

Since w — P2/P1, the limit w 00 corresponds to 
Pi = with p2 remaining nonzero. This means only chip- 
ping moves are allowed in the mass model. As mentioned 
in Section II, in this limit the model has some resem- 
blance to models of vacancy cluster formation |0 , with 
mass in our model representing the number of vacancies 
in a cluster. Chipping off from clusters with m > 2 cor- 
responds to vacancy detachment, whereas the chipping 
move for m = 1 is tantamount to either diffusive hopping 
(if the neighbouring site is empty), or to attachment to 
a cluster (if the neighbouring site has a mass on it). 

In the LG version, the 'only-chipping' model corre- 
sponds to just the exclusion process fl^,]!^ either sym- 
metric (SEP) (corresponding to SCMAM) or asymmet- 
ric (ASEP) (corresponding to ASCMAM). Both the SEP 
and ASEP have been studied in great detail and sev- 
eral exact results are known jl^JT^. In particular it is 
known that in the steady states of both SEP and ASEP, 
all configurations are equally likely [ p^|JT^ . This means 
that in the thermodynamic limit, the steady state has a 
simple product measure: the probability P{m) of having 
exactly m holes following a particle is simply given by, 
P(m) = p'{l - p'Y^. Using p' = 1/(1 + p), we thus ob- 
tain an exact result for P{m) of the mass model in the 
w 00 limit, 



P(m) 



{i + py 



(9) 



which clearly demonstrates the exponential decay of 
P{m) for large m. 

We note in passing that Eq.(^ also describes the dis- 
tribution of masses in any dimension in an 'only chipping' 
model in which the rate of chipping at site i is propor- 
tional to the mass nii on that site. In this case, evidently 
every unit of mass performs a simple random walk, and 
the problem is tantamount to that of M independent 
random walkers on the lattice. Standard methods of sta- 
tistical mechanics can then be used to describe the steady 
state. The grand partition function is then 



Z= (l + z + z^...) 



N 



1 



{1-z) 



N 



(10) 



where z is the fugacity, and the probability of finding m 
particles on a site is 



P{m) = (1 - z)z''' 



(11) 



Eliminating z in favour of p = z/(l — z), we see that 
Eq.(|ll]) reduces to Eq.(|). 



B. No chipping: w = 

We now consider the other limit w = where there 
is no chipping and the masses only diffuse as a whole 



and aggregate with each other. In the LG language, this 
would mean that particles undergo only long range hop- 
ping to the farthest available hole (without crossing the 
next particle). 

In order to find the exact steady state in this limit, we 
first consider the mass model on a finite ring i? of TV sites 
and total mass M — pN. As time progresses, masses dif- 
fuse and coagulate with each other and the number of 
empty sites decreases. Since the diffusive motion con- 
fined in a finite region of space is ergodic, eventually all 
the masses will coagulate with each other and there will 
be precisely one single big conglomerate with mass M 
which will move on the ring. Thus in the steady state, 
one gets exactly, P{m) = 5{m — pN). In the LG lan- 
guage, this would mean that on a finite lattice, the par- 
ticles and holes will become completely phase separated 
in the steady state. 

It is also useful to study the approach to this phase 
separated steady state. We first note that on an infinite 
1-d lattice, the time dependent single site mass distribu- 
tion function, P(rn, t) can be exactly solved in the w — 
limit ||2^. It was shown exactly in Ref. that in the 
scaling limit, m — s- 00, t ^ 00 but keeping mj \Jt fixed, 
the function P{m,t) ^ t^^^'^S{m/c\/t) where the con- 
stant c depends on the initial condition but the scaling 
function S{x) is universal and is given by 



Six) 



■ exp [- 



TTX 



(12) 



In the LG model, this would mean that as time pro- 
gresses, the system undergoes phase separation and 
breaks into domains of particles and holes. The average 
linear size of these domains grows with time as l(t) ^ t^/^ 
at late times. At this point it is useful also to note that in 
usual models of coarsening with locally conserved dynam- 
ics, the domain size grows as l{t) ~ t^^^ |2^ as opposed 
to t^^^ here. This, however, is not entirely surprising 
since in our model, even though the particle number is 
conserved globally, the long range hopping effectively re- 
duces this to a locally nonconserved model with growth 
law t^/'^. Similar behaviour was noticed earlier in other 
models of coarsening with globally conserved dynamics 
psf . In an infinite system the domain size keeps growing 
indefinitely as l{t) ^ t^^^. However in a finite lattice of 
sites L, when l{t) ^ L after a time t ^ L^, the domains 
stop growing and the system eventually breaks up into 
two domains only, one of particles and the other of holes. 
We note that all the conclusions made in this subsection 
are equally valid for both symmetric and asymmetric ver- 
sions of the model. 

To summarize this section, we find that for fixed p in 
one dimension, the two extreme limits w — > and w — s- cx) 
are exactly solvable for their steady states. In the limit 
w ^ 0, the system has a phase separated steady state 
in the LG model which in the mass model corresponds 
to having a single massive aggregate. In the other limit 
?z; 00, the LG corresponds to simple exclusion process 
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with product measure steady state, which corresponds 
to an exponential mass distribution in the mass model. 
Thus there is a competition between the long range hop- 
ping that tends to create phase separation and the short 
range hopping that tends to mix the particles and holes 
to produce a product measure steady state. As w is in- 
creased from for fixed p, the massive aggregate coexists 
with power law distributed smaller masses (or hole clus- 
ters in the LG language) up to some critical value Wc{p)- 
For w > Wc{p), the massive aggregate disappears and the 
cluster size distribution of holes becomes exponential. 

The phase separated state found in the lattice model 
is reminiscent of the state found in lattice gas models 
of traffic with quenched disorder, either in the distribu- 
tion of maximum velocities of different cars or at 
different locations on the road ||l^ , [l5|l . Our model, by 
contrast, has no disorder and is translationally invariant 
to start with, although this symmetry is broken in the 
phase-separated state. 



VI. DYNAMICS IN ONE DIMENSION: 
MAPPING TO INTERFACE MODELS IN (1 
DIMENSION 



In the previous sections we have studied the static 
properties of the model in the steady state. However it 
is also important to study the dynamics of the model in 
the steady state. The universal features of the dynamics 
is usually best captured by time dependent correlation 
functions in the steady state. In this section we how- 
ever take a slightly different route. Instead of studying 
the time-dependent correlation functions directly in the 
mass or the equivalent latice gas model, we first map the 
LG model onto an interface model and then study the 
time dependent properties of the width of the fluctuat- 
ing interfaces. The advantage of this route is that not 
only does it capture the essential universal features of 
the dynamics, but it also makes contact with other well 
studied interface models in certain limiting cases. 

There is a standard way p7[ | to map a LG configuration 
in one dimension to that of an interface configuration on 
a 1-d substrate. One defines a new set of variables {Si\ 
such that S'i = 1 if the i-th site is occupied by a particle 
and Si — —1 ii it is empty. Then the interface height hi 
at site i is defined as, hi — ^j- Thus the overall 

tilt of the interface, tan 6* — {h^ — hi)/ L = 2p' — 1 is set 
by the particle density p' in the LG model. 

The different dynamical moves in the LG model can 
be translated in a one to one fashion to corresponding 
moves of the interface. For example, if a particle at site 
i jumps one unit to the right, it corresponds to the de- 
crease of height hi by 2 units, hi ^ hi — 2. Similarly, if 
the particle at site (i -I- 1) jumps one unit to the left, the 
height hi increases by 2 units, hi — > hi+2. Thus the near- 
est neighbour hopping of particles to the left or right in 
the LG model corresponds respectively to deposition and 



evaporation moves in the interface model. Similarly the 
long range hopping of a particle to the farthest available 
hole before the next particle would translate to nonlocal 
moves in the interface as shown in Fig. 3. Once again, 
the ratio of the rates of evaporation-deposition moves to 
that of the nonlocal moves is given by the parameter w. 

A natural measure of the fluctuations of the interface 
is its width deflned for a flnite system of size L as. 



W{L,t) 



\ 



(13) 



where h = J2i=i ^i/L. The width W{L^t) is expected to 

(14) 



have a scaling form ||1J^ 

W{L,t) - L^F{t/L'), 



in the scaling limit with large L, large t but keeping 
t/L^ finite. The exponents x ^-nd ^ respectively the 
roughness and the dynamical exponents and the scaling 
function F{x) is universal with the aymptotic behaviour: 
F{x) 0(1) as a; ^ oo and F{x) ^ as x ^ where 
(3 — Xz. The exponents x ^^'^ z characterize the univer- 
sality classes of the interface models. 

In the following we keep the particle density p' fixed 
and investigate the width of the corresponding interface 
model W{L,t) by varying the parameter w. How do the 
exponents x ^-nd z that characterize the universal be- 
haviour of interface fluctuations change as one varies w 
from to oo? We study this question for both the sym- 
metric and the asymmetric models and denote the re- 
spective exponents by (Xw^'^K z^^'^'^) and {Xw'''^'''\ Zw^"''^^)- 
In the following three subsections, we study the width 
W{L,t) in both models for three special values of w, 
namely w oo, w = and w = Wc{p')- In the flrst 
two cases we find analytical results whereas in the last 
case, we present numerical results only. 



A. Only chipping: w ^ oo 

Wc first consider the symmetric mass model SCMAM 
whose lattice gas equivalent corresponds to the SEP in 
the w ^ oo limit as noted in section V. This means that 
in the corresponding interface model, nonlocal moves are 
absent and only the allowed moves are evaporation and 
deposition subject to certain local constraints. An im- 
portant property of this symmetric model is that the av- 
erage velocity of the interface is zero as the evaporation 
and deposition occurs with equal probability. It is well 
known [T7| that the continuum version of this discrete 
interface model corresponding to SEP is well described 
by the Edwards- Wilkinson equation ||l8|,|l9|. This linear 
evolution equation can be easily solved and one gets the 
exact exponents ]l8|,p^. 



(s) i. 

' " 2- 



(^) - 2 



(15) 
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The LG equivalent of the asymmetric mass model AS- 
CMAM in the w; ^ oo limit is the ASEP. In the cor- 
responding interface model the allowed moves are either 
only deposition or only evaporation (but not both) de- 
pending upon the direction of particle motion in ASEP. 
Thus the interface has a nonzero average velocity and 
its continuum version is known |17| to be described by 
the nonlinear KPZ equation ]20|] . In one dimension, the 
KPZ equation can be solved exactly and the exponents 
are known to be [E0| , 



Xc 



(as) _ 



1 

2' 



(as) 



(16) 



Thus in the limit w —^ oo, our model reduces to two 
known interface models for symmetric and asymmetric 
cases respectively and the corresponding exponents are 
obtained exactly. 



B. No chipping: w — 

We have mentioned in section V that in the w — > 
limit, as time progresses the LG phase separates into 
domains of particles and holes. In the equivalent spin 
representation where a particle is represented by an up 
spin. Si = 1 and a hole by a down spin, Si = —1, this 
then represents a spin model coarsening with time with 
average domain size growing as, l{t) ^ t^^^. In a finite 
system of length L, eventually the system breaks up into 
two domains of opposite signs. In the interface represen- 
tation, this would mean that the system would develop 
a single mound. Mound formation has also been studied 
recently in other interface models ||2^ . 

We note from Eq. |l3| and the definition, hi — 
that the calculation of the width W{L,t) requires the 
expression for the equal time spin correlation function, 
{S{i,t)S{j,t)) in the equivalent spin model. From the 
general theory of coarsening it is known ||2^ that for 
l{t) << L, the equal time spin correlation function sat- 
isfies the scaling behaviour, {S{0,t)S{r,t)) ~ G{r/l{t)). 
We have assumed that the system size L is large so that 
translational invariance holds in the bulk of the system. 
Using this scaling form in Eq. ^ and the result, l{t) ~ 
i^/^, a simple power counting gives, W{L, t) ~ LF{t/L'^). 
Thus for w = 0, we have the exponents. 



Xo = 1; 



zq = 2, 



(17) 



for both the symmetric and asymmetric models. The 
roughness exponent Xo = 1 can be easily understood 
from the fact that in the steady state the system de- 
velops a single mound whose maximum height is of 0{L) 
and hence the width of the fluctuations is also of order 
0(L). 



C. Critical Point: w = Wc{p) 

In this subsection, we keep the density p' ~ 1/(1 -|-p) 
fixed and tune w to its critical value, w = Wc{p) and cal- 
culate the width W{L,t) of the interface. Unfortunately 
we are unable to obtain any analytical result for this crit- 
ical case and will only present numerical estimates. 

We first consider the symmetric model. In this case 
we fix p = 1. In the LG language this means particle 
density, p' = 1/2. We first estimate the critical point, 
'Wc{p = 1) ~ 3.35 by simulating the equivalent mass 
model. Next we fix the value of p = 1 and w = 3.35 
in the interface model and measure the width W{L^t) 
for different lattice sizes, L = 16, 32, 64, 128 and 256. In 
order to verify the scaling form, W{L, t) ~ L^" F{t/ L^"), 
we plot in Fig. 4, W/L^" as a function of t/L^" for dif- 
ferent L. The best collapse of data is obtained with the 
choice, Xc^'^ ~ 0.67 and Zc^"^ « 2.1. 

The asymmetric case however is much more compli- 
cated due to the possible existence of logarithmic factors 
in one dimension and perhaps also due to other strong 
corrections to scaling. Proceeding as in the symmetric 
case, we first estimate the critical point, wdp = 1) ~ 0.77 
and the measure the width W{L,t). In this case we did 
not find a good data collapse using the canonical scaling 
form, W{L,t) ~ L^F{t/ L'^). Instead we tried collapsing 
the data assuming, W{L,t) WstF{t/ L^) where Wst is 
the steady state width W{L, oo). This gives a somewhat 
better convergence to collapse as L increases as shown 
in Fig. 5. with the choice z « 1.67. In the inset, we 
plot Wst vs. L. It is difficult to estimate the asymptotic 
growth Wst ^ (possibly with a logarithmic correc- 
tions) from the available data. A naive linear fit to the 
log-log plot of Wst vs. L gives an estimate of the slope, 
X ~ 0.68. Thus given the available data, we find the 
approximate estimates for the exponents, Xc-"'"'^ ~ 0-68 
and Zc*-"''^ ~ 1.67 at the critical point of the asymmet- 
ric model. However these are just approximate estimates 
and one needs larger scale simulations to determine the 
exponents more accurately for the asymmetric model. 



D. Flows 

We have also studied the width T/F(L,t) numerically 
for other values of w. This includes the numerical ver- 
ification of analytical predictions for the exponents in 
the limit of small and large w. We do not present here 
all these details but summarize the main picture that 
emerges from these studies by means of the schematic 
fiow diagram shown in Fig. 6. 

We find three different sets of exponents (x, z) that 
characterize the behaviour in three regions on the w axis 
for fixed p: subcritical when w < Wc{p), critical when 
w = Wc{p) and supercritical when w > Wc{p)- The sub- 
critical regime is controlled by the aggregation fixed point 
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(denoted AGG in Fig. 6) at w = 0, i.e., the phase sep- 
aration fixed point with x = 1 and z = 2 for both sym- 
metric and asymmetric models. The supercritical regime 
is controlled by the fixed point at w — > oo. For the sym- 
metric case, this is the Hammersley-Edwards- Wilkinson 
fixed point HEW (x — 1/2, z — 2), whereas in the 
asymmetric case this is the KPZ fixed point (x — 1/2, 
z = 3/2). The fixed points SG and ASG in Fig. 6 cor- 
respond to criticality in the symmetric and asymmetric 
conserved mass models respectively; these are new un- 
stable fixed points with exponents, (x ~ 0.67, z « 2.1) 
for the symmetric case and (x ~ 0.68, z « 1.67) for the 
asymmetric case respectively. 



VII. GENERALIZATION TO ARBITRARY 
FRAGMENTATION KERNEL: RELATION TO 
OTHER MODELS 

In the mass model discussed so far we have considered 
only two possible moves, namely "chipping" of a single 
unit of mass to a neighbouring site with rate w or hop- 
ping of the mass as a whole to a neighbouring site with 
rate 1. However in a more general setting, k units of mass 
can break off a mass m and hop to a neighbouring site 
with rate p{k\m) where k < m. In the equivalent lattice 
gas version in one dimension, this would mean the hop- 
ping of a particle with rate p{k\m) to the fc-th hole to 
the right or left without crossing the next particle which 
is located at a distance m + 1. In the asymmetric ver- 
sion, the particles jump only along one direction as usual 
(either left or right). In the model studied so far. 



p{k\m) = wSk,i + Sk,m, 



(18) 



where Sij is the Kronecker delta function. With this 
mass fragmentation kernel, we have seen that the in 
the steady state the system undergoes a nonequilibrium 
phase transition as the parameter w that controls the 
relative strength of the two delta peaks is varied. 

The question naturally arises as to whether this 
phase transition exists for arbitrary fragmentation ker- 
nel p{k\m). In general, it is hard to find the steady state 
analytically for arbitrary p(fc|m). However, for some spe- 
cial choices of the kernels, it is possible to obtain exact 
steady states via mapping to some other solvable models 
of statistical mechanics. We list below a few of them. 

In the context of traffic models, Klauck and Schad- 



schneider recently studied |26| an asymmetric exclusion 
process in one dimension where a particle can jump ei- 
ther to the neighbouring hole to the right with rate pi or 
to the second hole to the right with rate p2 ■ The corre- 
sponding jump kernel can be written as 



p{k\m) =piSkA +P26k,2 



(19) 



have simple product measure for all pi and p2- Using 
this result, it is easy to show that in the corresponding 
mass model with mass density p, the steady state single 
site mass distribution P{m) (same as the probability of 
having a hole cluster of size m in the LG model) is sim- 
ply given by, P{m) = p/{l + p)'"+^ and therefore decays 
exponentially for large m for arbitrary pi and p2. 

Recently Rajesh and Dhar studied |2^ an anisotropic 
directed percolation model in 3 dimensions. Their model 
can be reduced to an asymmetric hard core LG model 
with the following jump kernel. 



p{k\m) = /-'^'=«(1 -p)' 



i—k 



(20) 



where < p < 1. By mapping it to the five vertex 
model, the steady state of this model was shown ex- 
actly to have a simple product measure for all p. This 
then immediately gives for the mass model, once again, 
P(m) = p/(l + /9)™+i. 

Another exact result can be derived for the following 
asymmetric mass model. Instead of discrete mass, we 
now consider continuum masses at each site. The mass 
rrii at each site i evolves in discrete time according to the 
following stochastic equation, 

mi{t + 1) = qi-i,imi^i{t) + (1 - qi^i+i)mi{t) (21) 

where the random variable qi-i^i represents the fraction 
of mass that breaks off from site (z — 1) and moves to 
site i. We assume that each of these fractions are inde- 
pendent and identically distributed in the interval [0, 1] 
with some distribution function ri{q). We show in the 
Appendix that for the special case of uniform distribu- 
tion, i.e., 'q{q) = 1 for all q in [0, 1], the exact steady state 
distribution P{m) is given by, 



4m 

Pim) = -^e 
P 



-Ira I p 



(22) 



By generalizing the matrix product ansatz used for ASEP 
pO|, the steady state of this model was shown pffl to 



where p = mP{m)dm is the conserved mass density 
fixed by the initial condition. In the LG language, this 
would mean that if rui is the distance between i-th and 
(z -|- l)-th particle, then the i-ih. particle can jump to 
any distance between k and k -\- dk io the right (without 
crossing the next particle on the right) with uniform rate 
p{k\'mi)dk = -^dk. In order to verify the exact formula 
for P{m) in Eq. ^ we performed numerical simulation 
in 1-d. In Fig. 7, we show for p = 1, the perfect agree- 
ment between the theoretical cumulative mass distribu- 
tion, F{m) = /" P{m')dm' = 1 - e'^"" - 2me~^"^ and 
the numerically obtained points from the simulation. We 
note that a similar result was recently derived pjl in the 
context of the generalized Hammersley process |l32| . 

Thus in all these cases studied in this section, we find 
exponential decay of the mass distribution for large mass 
in the steady state. It would therefore seem that only 
for the special kernel, p{k\m) = w6k,i + Sk.m, is there 
a nontrivial phase transition from a phase where P{m) 
decays exponentially to another where it has a power law 
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decay in addition to an infinite aggregate. The presence 
of two delta function peaks in p(k\m) seem to be cru- 
cially responsible for this phase transition. We believe 
that the phase transition will still persist if one allows 
for a nonzero width to the delta peaks at the two ends of 
the kernel near k = and k = m but making sure that 
the widths remain finite even as m ^ oo. This however 
needs further studies to be confirmed. 



VIII. CONCLUSION 

In this paper we have studied a lattice model of ag- 
gregation and dissociation where a mass from a site can 
either move as a whole to a neighbouring site with rate 
Pi or can chip off a unit mass to a neighbour with rate 
P2- The hopped mass then aggregates instantaneously 
with the mass that is already present at the neighbour. 
The ratio of the rates, w = P2IP1 and the conserved mass 
density p are the only two parameters of the model. The 
steady state of the model undergoes a phase transition 
as the parameters (p, w) are varied. In the (p, w) plane 
there is a critical line that separates two phases: (i) the 
'Exponential' phase where the single site mass distribu- 
tion P(m) decays exponentially for large m and (ii) the 
'Aggregate' phase where P(m) has a power law decay 
in addition to a delta function peak at m = 00 signify- 
ing the presence of an infinite aggregate. On the critical 
line, the aggregate vanishes but P{m) ^ m^'^ retains the 
same power law tail. We have also studied how the uni- 
versality class of this dynamical phase transition changes 
on applying a bias in a particular direction of the mass 
transport. 

We have solved the model exactly within the mean field 
theory and presented numerical results for one dimension. 
Besides, by exploiting a mapping to a lattice gas model 
in 1-d, we have obtained the steady state distribution 
P(m) exactly for small and large w. We have further 
mapped the lattice gas model to an interface model in 
(1-1-1) dimension and studied the width of the interface 
that characterizes its fluctuations. We have calculated 
the roughness exponent x and the dynamical exponent z 
analytically for small and large w both with and without 
bias. We have computed these exponents numerically at 
the critical point and shown that the two critical points 
(with and without bias) represent two new fixed points 
of interface dynamics in (1 -t- 1) dimension. 

We have also generalized our model to include arbi- 
trary fragmentation kernels and obtained a few exact re- 
sults for special choices of these kernels via mappings to 
other solvable statistical mechanics models studied in the 
recent past. 

Our model obviously has some shortcomings as a 
model of realistic aggregation and fragmentation phe- 
nomena. For example, we have assumed that the rate 
of hopping of a mass as a whole is independent of the 
mass. In a realistic setting this is perhaps not a good ap- 



proximation. For realistic modelling of a specific system, 
one needs to put in these details. This however was not 
attempted in this paper, whose aim is to understand the 
mechanism of the dynamical phase transition induced by 
the basic microscopic processes of diffusion, aggregation 
and fragmentation within a simple setting. 

However there remain many open questions even for 
this simple model. For example, in one dimension we 
were not able to compute analytically the various expo- 
nents at the critical point. Given that there has been 
recent progress in determining the exact steady states of 
a class of exclusion processes in one dimension via using 
matrix product ansatz | pO[ | and coordinate Bethe ansatz 
p3[ , it may be possible to obtain the exact steady states 
of our model in 1-d. 

Another important open question that needs to be 
studied for both the symmetric and asymmetric models: 
What is the upper critical dimension dc of these mod- 
els beyond which the mean field exponents will be exact. 
From general analogy to other diffusion limited models 
studied earlier |3j], one expects that in the symmetric 
case dc — 2. This expectation is supported by the nu- 
merical fact that the the exponent r^. w 2.33 in 1-d is 
already close to its mean field value 5/2. This conjecture 
however requires a proof. 

The question of upper critical dimension is however 
puzzling for the asymmetric model where the numerical 
value of Tas ~ 2.05 (perhaps Ts — 2 with logarithmic 
corrections) is quite far from the mean field value 5/2. 
This indicates that the dc (if there exists one) for the 
asymmetric case might even be bigger than 2 which how- 
ever is quite contrary to the naive expectation that the 
dc of directed models is usually lower than that of undi- 
rected models. A numerical study in 2 dimensions of the 
asymmetric model may shed some light on this puzzle. 

We thank D. Dhar and Rajesh R. for useful discussions. 



APPENDIX: EXACT RESULT FOR UNIFORM 
FRAGMENTATION KERNEL 

In this appendix we solve exactly the steady state sin- 
gle site mass distribution function P(m) for the model 
where the mass rrii at every site i evolves in discrete time 
according to the stochastic equation, 

mj(t + 1) = qi-i,im.i^i{t) + (1 - qi^i+i)mi{t) (23) 

where the fractions g^.j are independent and identically 
distributed random variables in [0, 1] with distribution 
function ri(q). We note the formal similarity between 
Eq. 23 and the force balance equation in the g-model of 



Coppersmith et. al. 



(24) 



where W{i, D) represents the net stress supported by a 
glass bead at a depth Z? in a cylinder and qji is the frac- 
tion of the stress transported from particle j at layer D 
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to a particle i at layer {D + 1). The only difference be- 
tween Eq. [2^ and Eq. |^ is in the additional constant 
term 1 in Eq. |2j that is absent in Eq. |2^ . Nevertheless 
the same line of argument as in ref. |28|] leads to an exact 
solution in our case also as we outline below. 

We first consider the mean field theory where we ne- 
glect correlations between masses. For simplicity, we use 
the notations, = qi, l — qi^t+i = 92, = mi and 

rrii = 7712 ■ Mean field approximation then leads to to the 
following recursive equation for the mass distribution. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 

[7] 



P{m,t+ 1) 



1 /•! poo pQO 

/ dqidq2i]{qi)7]{l - q2) / / dr7iid77i2 [8] 
Jo Jo Jo 



X P{mi,t)P{77l2,t)S{777 - 77liqi - r7l2q2), (25) [9] 



where we have used Eq. |2^. In the limit t ^ 00 and for 
the uniform distribution 7]{q) = 1, the Laplace transform 
P{s) of the distribution P{m) satisfies the equation 



P{s) 



dqP{sq)\ 



(26) 



Defining V{s) — y P(s) and u — qs, we get from Eq. |2q , 



sV{s) 



duV'^{u). 



Differentiation with respect to s yields, 
as 

which can be integrated to give, 

1 



Vis) = 



1 - Cs 



(27) 



(28) 



(29) 



The constant C is determined from the mass conservation 
equation, mP{77i)d77i — —dP{s)/ds\s^o — p where p 
is the conserved mass density. Thus, C — dV/ds\s=o = 



-pl2. Hence we get, P(s) = 1/(1 
Laplace transform, 



P(to) 



4™ 



-2m/ p 



\s) and by inverse 



(30) 



This is the mean field result for P(m). However in ref. 
p8[ , it was proved that for uniform distribution 77(9) = 1, 
the mean field stress distribution, where the stresses sat- 
isfy Eq. |2^, is exact. The explicit algebraic proof in ref. 
psf proceeded via constructing exact recursion relations 
for the joint probability distribution of the weights in 
row {D + 1) in terms of those for row D, and showing 
that the mean field factorization of these joint distribu- 
tions are invariant under this recursion. The same line of 
proof can be adapted to show that the mean field result 
Eq. ^ is also exact for our problem. We do not giv e 
details of the proof here as they can be found in ref. [E8| . 
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FIG. 1. log-log plot of P{m) vs. m for SCMAM for two choices of the parameters: (u; = 1.0, p = 0.2) and (w = 1.0, p = 3.0) 
on a periodic lattice of size A'^ = 1024. Inset: Phase diagram. The solid line and the points indicate the phase boundary within 
mean field theory and 1-d simulation respectively. 




FIG. 2. log-log plot of P{m) vs. m for ASCMAM for two choices of the parameters: {w = 1.0, p = 0.2) and {w = 1.0, p = 10.0) 
on a periodic lattice of size N = 1024. Inset: Phase diagram. The solid line and the points indicate the phase boundary within 
mean field theory and 1-d simulation respectively. 
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FIG. 3. The constructions of the equivalent lattice gas model and interface model are illustrated for a particular configuration 
of the conserved mass model. The shaded blocks indicate the masses that would move as a result of "diffusion and aggregation" 
and "chipping". 



14 



0.25 



0.2 



w L 



0.15 





L=16 + 




L=32 X 




L=64 I 








L=266 ■ 










's 


- 


1 




1 




B 





0.05 



0.15 



tL 



FIG. 4. The scaled width W/L^ is plotted against the scaled time f/L^ for lattice sizes L = 16, 32, 64, 128 and 256 at the 
critical point (p = 1, Wc ~ 3.35) of the interface corresponding to SCMAM. The best data collapse is obtained with the choice 
of exponent values, x 0.67 and z m 2.1. 
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FIG. 5. The scaled width W/Wst is plotted against the scaled time t/L^ for lattice sizes L — 16, 32, 64, 128 and 256 at 
the critical point (p = l,uic « 0.77) of the interface corresponding to ASCMAM. The best convergence to data collapse as L 
increases is obtained with the choice z « 1.67. Inset shows Wst plotted against L. 
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FIG. 6. Schematic depiction of fixed points and associated flows 
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FIG. 7. The theoretical cumulative mass distribution function, F(ni) — 1 — e ^"^ — 2me ^"^ (continuous line) plotted against 
m along with the numerically obtained points from the simulation of the model in one dimension for mass density p = 1. 
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